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ABSTRACT 

Torsion oscillations of the neutron star crust are Landau damped by the Alfven 
continuum in the bulk. For strong magnetic fields (in magnetars), undamped Alfven 
eigenmodes appear. 



1. Introduction and Conclusion 

It is thought that Israel et al (2005) have detected torsion oscillations of the neutron star crust. 
Here we show that crustal modes are strongly affected by the Alfven continuum in the bulk. For 
week magnetic fields, there is Landau damping of crustal eigenmodes. At stronger magnetar-like 
magnetic fields undamped Alfven eigenmodes appear. In this regime, the Alfven waves in the bulk 
control the torsion oscillations of the crust. 

It might be that Israel et al (2005) have actually measured the Alfven speed in the bulk 
ca, rather than the torsion speed in the crust q. If so, the measured v ~ 20Hz frequency gives 
c A « 2uR ps 4 x 10 7 cm/s (assuming that v is the fundamental Alfven eigenmode). 

We calculate the Landau damping of the crustal torsion modes and Alfven eigenmode (un- 
damped) frequencies. The thin-crust uniform-field constant-density model is used, and we only 
consider the antisymmetric axisymmetric modes. This is obviously not good enough for serious 
neutron star seismology. One must include: (i) the non-axisymmetric modes, (ii) the toroidal mag- 
netic field, (iii) the realistic density profile. But the qualitative results - Landau damping and 
emergence of undamped Alfven eigenmodes - should be robust. 

This paper is a formal mathematical development of the ideas of Levin (2007). 



2. Results 

We use the following units and notations: 

• R = 1 is the stellar radius 

• ca = 1 is the Alfven speed in the bulk 
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• 3M is the crust to bulk mass ratio. 

• ct is the torsion sound speed in the crust (in units of ca)- 

The calculated frequencies u> (in units of ca/R) of the antisymmetric axisymmetric torsion 
modes (I = 2,4,6, m = 0) are given in the table for various mass ratios M and speed ratios q. 
The eigenmodes (and cuts, see §3 ) are also shown in the figure for M = 0.1, c t = 2. 

3. The eigenmode equation and the dispersion relation 

Toroidal displacements of the crust launch Alfven waves into the bulk. For an axisymmetric 
antisymmetric torsion mode, the toroidal displacement in the bulk is oc e~ lujt sin(u;z), where z is 
the Cartesian coordinate along the magnetic field, with the equatorial plane z = 0. Including the 
corresponding force into the thin-crust eigenmode equation, we get 

c? ((1 - , 2 )/" - 4xfj + u, 2 / = "s^"*) ,. (i) 

Here x = cos 6, 6 is the spherical angle with 9 = along z, fsinO is the toroidal displacement of 
the crust, prime is the x-differentiation. The detailed derivation is given in the Appendix. 

For the antisymmetric mode, the boundary condition at x = is / = 0. The boundary 
condition at x = 1 is given by the equation itself: —4c 2 /' + J 1 f = i ^jf^-f- 

To obtain the dispersion relation in the form F(uj) = 0, choose some /(l) and /'(l) satisfying 
the x = 1 boundary condition. Then integrate (1) form x = 1 to x = and put F(cj) = f(0). 

According to the Landau rule, one calculates F{oS) in the upper half-plane of complex to. The 
eigenmodes are given by the zeros of F(cj) on the real axis (undamped modes) and in the lower 
half plane (Landau damped modes). The value of F in the lower half-plane should be obtained by 
the analytic continuation from the upper half-plane. 

Table 1: Eigenfrequencies in units of ca/R- 



M 


= 0.1, c t 


- 1 


2.10, 2.74, 3.03, 3.12, 5.05-0.001i, 5.75-0.0021, ... 


M 


= 0.1, c t 


= 3 


3.06, 6.18-O.lli, 9.35-0.041, ... 


M -- 


= 0.1, c t 


= 10 


21.7-4.01i, 46.3-6.19i, ... 


M = 


~- 0.1, c t = 


= 100 


200-3.14i, 424-3.M, ... 


M = 


-- 0.01, c t 


-= 10 


3.11, 6.26-0.005i, 9.41-0.006i, ... 


M = 


-- 0.01, c t 


= 30 


103-76i, 198-144i, ... 


M = 


0.01, c t 


- 100 


219-551, 490-102i, ... 
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The analytic continuation of F to the lower half-plane can be accomplished by complexifying 
x and integrating (1) from x = 1 to x = along an arbitrary contour in the upper half-plane of x. 
The shape of the x-contour determines the shape of the cuts in the lower half-plane of uj. We chose 
the semicircle \x — 0.5| = 0.5. This gives the cuts in the lower half-plane of uj running straight down 
from uj = irk. 

The dispersion function F was calculated numerically. The eigenmodes, both the Landau- 
damped crustal modes and the undamped Alfven eigenmodes are given in the table. The real 
entries of the table (what we call Alfven eigenmodes), were confirmed by a straightforward real- 
numbers integration of (1). 

Due to the cuts, the late-time asymptotic of the crustal motion will also have algebraically 
damped modes with frequencies uj = nk. The late time asymptotic is determined solely by the 
location of the tips of the cuts and is not affected by the choice of the x-contour. The zeros of F 
are also independent of the x-contour. 

For generic parameters M and q, numerical integration seems to be the only way. But there 
are limiting cases which can be treated analytically. These may serve to confirm that equation (1) 
actually makes sensible predictions and also to check the numerical results: 

• q « 1, M » 1 gives undamped crustal modes uj 2 = (l 2 + l — 2)c£ + I = 2,4,6, .... This case 
is obvious - magnetic field and the fluid bulk just follow the slow and heavy crust, providing an 
additional elasticity and therefore increasing the frequency. This case is unphysical, because 
the crust is actually lighter than the bulk. 

• ct 3> 1, Met ~> 1 gives Landau damped crustal modes uj = Vl 2 + 1 — 2 ct — ijfc, I = 2,4, 6, 

ai are calculable dimensionless numbers (of order unity for low I). This case is physical, it 
occurs for small enough magnetic fields. 

The calculation is as follows. By the Euler's formula, 



We solved it using the first order perturbation theory in 1/(Mq). For / = 2, one calculates 
0-2 = 5/16 = 0.313 - in agreement with the M = 0.1, q = 100 entry of the table. 

I thank Yuri Levin for showing me the problem and for useful discussions. 
This work was supported by the David and Lucile Packard foundation. 





Now (1) becomes 




UJX 



(3) 
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A. Torsion oscillations of the elastic crust with the magnetized bulk 

Consider a neutron star with a thin crust. Choose the z axis along the uniform magnetic field 
B. The toroidal displacement of the fluid bulk £ is described by the Alfven wave equation 

dU = <? A d% (Al) 
where c\ = is the Alfven speed in the bulk of density p. 

The toroidal displacement of the crust ip is described by the torsion wave equation with the 
magnetic force from the bulk 

(A2) 

Here q is the torsion sound speed in the crust, R is the radius of the star, 6 is the polar angle with 
respect to z, T^ z is the toroidal-vertical component of the Maxwell stress tensor, a is the surface 
density of the crust. 

The Maxwell stress is given by T^ z = —3^, B^ = Bd z ^, which should be calculated at the 
boundary of the bulk, where £ = ip. We get: 

^ = 4 (%1> + ~ -^r-J + 2^ - ^ cos ed z ^\ z=cosd . (A3) 

R 2 \ sm9 sir 9 ) a 

As a check, one confirms that this system conserves energy and angular momentum: 

E = \j pdV ((d t 2 + c\(U) 2 ) + \J 'odA {{d^f + -|( s in^(-^)) 2 ) (A4) 

L = j pdV rsinO d t £ + J odA RsinO d t ip, (A5) 
where J dV and J dA are volume and surface integrals and r is the radius. 
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Fig. 1. — M = 0.1, ct = 2. Zeros and cuts of the dispersion function F(u>) in the complex uj plane. 



